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We present results for the phase diagram of an SU(JV) generalization of the Heisenberg antifer- 
romagnet on a bipartite three-dimensional anisotropic cubic (tetragonal) lattice as a function of TV 
and the lattice anisotropy 7. In the "isotropic" 7=1 cubic limit, we find a transition from Neel to 
valence bond solid (VBS) between TV = 9 and TV — 10. We follow the Neel-VBS transition to the 
limiting cases of 7 <C 1 (weakly coupled layers) and 7 ^> 1 (weakly coupled chains). Throughout the 
phase diagram we find a direct first-order transition from Neel at small- TV to VBS at large-TV. In 
the three-dimensional models studied here, we find no evidence for either an intervening spin-liquid 
"photon" phase or a continuous transition, even close to the limit 7 1 where the isolated layers 
undergo continuous Neel-VBS transitions. 



I. INTRODUCTION 

The theoretical and experimental study of the destruc- 
tion of magnetic order at T = has become one of 
the most popular topics of quantum condensed matter 
physics in recent decades.^ Such phase transitions arise 
in a general class of many body systems referred to as 
"quantum magnets," prominent examples being transi- 
tion metal oxide systems, heavy fermion materials, and 
ultra cold atomic gases in optical lattices. Given the 
great interest in this phenomena in experimental systems, 
an extensive body of theoretical work has been devoted 
to the study of the simplest microscopic models that dis- 
play these phase transitions, i.e., quantum spin systems. 
Model systems have been studied theoretically by a large 
variety of analytical and numerical methodsPHS 

The most popular spin models are anti-ferromagnets 
with SU(2) symmetry. Later generalizations to SU(TV) 
symmetry have also received a lot of attention, both in 
the large-TV limit and for finite values of TV. For bipar- 
tite lattice models, the simplest magnetic order is of the 
semi-classical kind, which we will term Neel order for 
all TV. On increasing quantum fluctuations, theoretical 
reasoning suggests that a natural non-magnetic state to 
consider in the case of many anti-ferromagnets is a trans- 
lational symmetry breaking - valence bond solid (VBS). 
A natural line of investigation that then arises is the the- 
oretical description of the Neel-VBS transition. Inter- 
estingly the nature of this transition is strongly affected 
by the dimensionality of the system. In this work we 
study the transition between Neel and valence bond solid 
(VBS) phases in antiferromagnets with SU(TV) symmetry 
in three spatial dimensions by unbiased quantum Monte 
Carlo (QMC) simulations. Before we turn to an exposi- 
tion of our results, we will briefly review past results of 
studies of the SU(TV) Neel-VBS transition in one and two 
dimensions. 

In one dimension the transition between "Neel" and 
VBS can be studied accurately by controlled theoretical 
methodsP We used the quotation marks because in one 
dimension there can be no true long-range order for the 
Neel phase, although true long-range order is possible at 
T = for the VBS. The transition between the power- 




FIG. 1. (color online) . Three dimensional lattice on which the 
model, Eq.[l] is defined. Left: A 2x2x2 piece of the cubic lat- 
tice showing nearest neighbor (Ji and 7J1) and next-nearest 
neighbor (J2 and 7J2) couplings. Bonds in the x — y plane 
have strength Ji , J2 . The parameter 7 is a multiplicative fac- 
tor for all bonds with a ^-component that connect neighboring 
x — y planes. Right: The pattern of red bonds depicts one 
of the degenerate columnar dimer patterns found in the VBS 
phase of the model given by Eq. [l] The dimer pattern is such 
that the Bragg peaks are located at (n, 0, 0) and (0, n, 0) for 
7 < 1 and at (0,0, 7r) for 7 > 1. At 7 = 1 all three Bragg 
reflections appear. 



law "Neel" and VBS states has been studied in detail 
in the anti-ferromagnetic SU(2) J\-J<i model.^ More re- 
cently, the same transition has been s tud ied in sign prob- 
lem free one-dimensional JQ models!^ The theoretical 
expectation of a Kosterlitz-Thouless type transition be- 
tween "Neel" and VBS states has been verified in numer- 
ical studies. 

In two dimensions the situation gets more interest- 
ing since true long range order is possible for the Neel 
state at T = 0. While conventional Landau theories 
would predict a first-order transition or an intervening 
phase, a field theoretic work made the prediction of a 
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direct continuous phase transition between conventional 
Neel ordered and VBS phases, called a deconhned crit- 
ical point There is now strong evidence that the de- 
confined scenario is realized in a series of models that 
allow the transition to be studied by QMC simulations 
as a function of N in SU(JV) antiferromagnetsP-^^ The 
SU(iV) studies give quantitative estimates for scaling di- 
mensions of both Neel and VBS order parameters for 
large values of N. By comparing these universal num- 
bers to the estimates from a 1/N expansion^ f the 
relevant thr ee-dim ensional classical non-compact CP^ -1 
field theoryE^SD w hich had previously been computed 
analyticallyp^flSl substantial credence was given to the 
identification of deconfined criticality in this system. 

In three dimensions the nature of the Neel- VBS phase 
transition has been studied only in SU(2) models^ with 
full cubic symmetry, where clear evidence for a first-order 
transition was reported. In this work we extend the pre- 
vious study^Slin two ways of interest: we study the Neel- 
VBS transition in systems that have SU(7V) symmetry 
with 2 < N < 10, and we study lattices with anisotropic 
couplings so that we can study the phase transition in 
the limit of weakly coupled two-dimensional (2D) layers. 
The motivation for these generalizations comes from our 
theoretical understanding of the Neel- VBS transition in 
two dimensions, where deconfined criticality arises due to 
the instability of three-dimensional (3D) compact gauge 
theories to confinement. In (3 + 1) dimensions, no in- 
stability exists for a compact gauge theory and thus one 
might expect to find an exotic gapless "photon" phase 
close to a transition between Neel and VBS. 

In this work, we endeavor to answer the following ques- 
tions: If one were to couple together 2D planes, each 
residing close to the deconfined quantum critical point, 
would the emergent gapless photon remain deconfined 
over some finite region of the phase space realizing the 
eagerly sought after photon phase in three dimensions? 
And, if there were no intervening phases between the Neel 
ordered and VBS phases, what is the nature of the tran- 
sition between the two phases in three dimensions? The 
search for a photon phase in microscopic models has been 
carried out in various contexts. Models in which there is 
evidence for the photon phase incl ude la rge- N versions 
of Sp(iV) magnets^ dimer modelsp^S] and U(l) sym- 
metric magnets pSM] 

In our study here, we will work on a bipartite lattice in 
which the spins transform as a fundamental representa- 
tion on one sub-latti ce and a conjugate to fundamental on 
the other sublattice ! 27 * 28 ! The specific lattice considered 
here may be thought of as an anisotropic cubic lattice 
(really a tetragonal lattice), where the coupling on the 
bonds in one direction (say the z direction) is different 
from the coupling on the bonds in the other two direc- 
tions [x and y directions). For SU(2) it is well known 
that independent of the anisotropy of the lattice, this 
model is Neel ordered. On the other hand, when N is 
very large the SU(iV) antiferromagnet maps to a quan- 
tum dimer model with only a kinetic term {22 This model 
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FIG. 2. (color online). Phase diagram of the Ji model, i.e., 
Eq. |TJ with J2 = 0, showing the transitions between Neel and 
VBS order as 7 is varied for values of N ranging from 2 to 
10. For N — 5, 6, 7, 8, 9, there are three segments for each N: 
the first (gray) corresponds to VBS order, the second (blue) 
corresponds to Neel order, and the third (gray) is again VBS 
order. For N = 2, the system is always magnetically ordered 
even in pure ID (the so-called Bethe-chain phase). Likewise, 
for N = 10, the system is always VBS ordered in three or less 
dimensions. Also note that for N = 3, 4, there are only two 
segments, Neel then VBS, as the system remains Neel ordered 
on the 2D side of the phase diagram since N — 5 is required 
in pure 2D to destroy the magnetic order. The horizontal 
axis is shown on a log2 scale. The valu e of N for which the 
system VBS orders in the 2D limhP^ and ID limilP 1 have 
been studied previously. 



is valence bond solid ordered at T = independent of the 
extent of anisotropyP^ From these simple considerations, 
it is clear that regardless of the anisotropy there must be 
a transition between Neel and VBS at some finite value 
of N. The central questions we have investigated here 
are whether there is a new intervening phase between 
the Neel and VBS phases and the nature of the transi- 
tions, i.e., whether they are first order or continuous as 
a function of N. 

In Sec.|nJ we define a 3 D mo del, by generalizing the 2D 
J\-Ji modelp^and in Sec. Ill we explain the methods by 
which we simulate our model. In Sec. |IV| we present the 



phase diagram inferred from our QMC simulations, and 
in Sec. [V] we study the nature of the phase transitions. 
Finally, in Sec. |VI| we summarize our results and present 
an outlook on future work. 



II. THE MODEL 

We construct our model from simple SU(iV) invariant 
operations. We first imagine our system of spins on some 
bipartite lattice is restricted to a Hilbert space wherein 
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the states on the sites of the A sublattice transform un- 
der the fundamental representation of SU(iV) while the 
states on the sites of the B sublattice transform under 
the conjugate to the fundamental represe ntation, a con- 
vention common in earlier works! 31 ! 32 ! 3 ' This results in 
a state of the form Y^ Q \a) A \a} B being an SU(iV) sin- 
glet. The first simple operation we include in our Hamil- 
tonian is the projector onto this singlet between 
two sites i and j on different sublattices. If one takes 
H = —Ji/Nj2(ij)Pij7 this reduces to the usual near- 
neighbor Heisenberg antiferromagnet (up to a constant) 
for N = 2. On the square lattice in two dimensions, it 
is already known that the Neel order gives way to VBS 
order for integer N > 5 (Refs. I3TI and I31?|): hence, the 
transition occurs between N = 4 and N = 5 if N is 
viewed as a continuous parameter. We expect that on 
the cubic lattice in three spatial dimensions, there will 
similarly exist some value of N (presumably larger than 
5) for which the Neel order fails to persist. We include a 
second term of another simple operation that permutes 
the spins on two sites of the same sublattice, Uij. This 
term is the SU( N) generalization of the Heisenberg ferro- 
magnet. With only this second term present, the system 
will trivially order ferromagnetically on each sublattice 
for any value of N. Turning on the first, antiferromag- 
netic term with small coupling will lock in the Neel or- 
der. Therefore, if we start with only the first term and 
some N sufficiently large such that the magnetic order is 
destroyed, we can presumably turn on and increase the 
coupling of the second term until Neel order is restored. 
In this way, we can continuously tune to the quantum 
transition point for each N on the VBS side of the tran- 
sition. See Fig. [I] for a depiction of the lattice and the 
relevant couplings. 

Finally, we would like to be able to vary the strength of 
the coupling between the 2D layers so we can interpolate 
from the results already known in one and two dimensions 
to the isotropic 3D limit. We accomplish this with the 
anisotropy parameter 7. Denoting the vector pointing 
from site i to j as fy, we have: 

E p,-4 E p - w 

rij-z=0 Vij-z^Q 

~§ E n «-r^ E n** 

((H)) ((H)) 

where (ij) denotes nearest neighbors and ((ij)) denotes 
next-nearest neighbors on the underlying cubic lattice. 
The limit 7 = recovers the 2D case, 7 — ► 00 gives 
a system of decoupled one-dimensional (ID) chains, and 
7 = 1 is the 3D isotropic case. It is worth mentioning here 
that the situation in ID is also well studied; assuming 
J2 = 0, only for N — 2 do we have the analog of Neel 
order, the so- called Bethe phase; for N > 3, there is again 
the VBS order2323! indicating a transition between N = 2 
and TV = 3. One of the goals of this paper is to determine 
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FIG. 3. (color online). Plots demonstrating the disappear- 
ance of magnetic order and onset of VBS order. The data 
shown here are for 7=1. (a) Spin stiffness, p„ (logio scale), 
vs. l/L for various values of N. (b) Square of VBS order 
parameter, xv, vs. l/L for various values of N. 

which integer values of N correspond to Neel and VBS 
order as we vary 7 from zero to infinity. The coupling 
J2 is introduced so that we can tune continuously across 
the Neel- VBS phase transition at fixed values of N and 
7- 



III. METHODS AND MEASUREMENTS 

All of our measurements were obtained exactly (to 
within statistical error) using the stochastic seri es ex- 
pansion QMC method with global loop updates! 4 ! 35 ! 36 ! 
The model presented in the previous section explicitly 
has all non-positive matrix elements leading to positive 
weights in the partition function and therefore no sign 
problem. Finite size scaling is performed on each quan- 
tity using a series of system sizes that varies depending 
on the anisotropy parameter. The system size is always 
taken to be L x L x L z where the ratio L/L z is adjusted 
to accommodate the anisotropy in a manner that will be 
explained below. 

Neel order is detected via the measurement of a finite 
spin stiffness p s , which is determined from the fluctua- 
tions of the spatial winding number of the world lines, 
Wf,. Specifically, PL^p s ^ = ((W^) 2 ), where \i = x,y,z 
and L x = L y = L. In the isotropic case, we report 
the average of all three components of the spin stiff- 
ness, p s = (p s , + Ps, z )/3 while on the 2D side 
of isotropy, we only average the x and y spin stiffnesses: 

Ps = (ps,x + Ps,y)/2. 
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VBS order is detected first by measuring the bond 
correlations corresponding to the predicted VBS pattern 
(in this case, columnar). The bond operator is simply 
Pij = P r< , rj (r) as defined in Sec. |llj The relevant bond 
correlator is thus: 



CV(r,r) 



N 2 



(P o ,x(0)P r , r+4 (r)> - <P o ,i(0)) 2 . (2) 



This does assume that the bonds will be oriented in the re- 
direction; in general, there is a degeneracy of the ground 
state that depends on the value of 7 (4 for < 7 < 1; 6 
for 7=1; and 2 for 1 < 7 < 00) and we measure the same 
correlator for y- and z-directed bonds. We average the 
x and y correlators for 7 on the 2D side of isotropy; we 
average all three directions at exactly 7 = 1; and we only 
look at the z correlator for 7 on the ID side of isotropy. 
Fourier transformation of this function and evaluation at 
uj = and k^, = np, (where fi = x, y, z is the relevant 
spatial coordinate) yields the Bragg peak corresponding 
to the VBS order. This quantity diverges as the quan- 
tum volume in the thermodynamic limit. If we divide by 
f3L 2 L z we get the square of the order parameter for the 
VBS phase, which we call xv- 



1 



Xv 



PL 2 L : 



rP 

/ drVCy(r,r)e J 
Jo 



(3) 



Another signature of VBS order is a finite energy gap 
for L — > 00. This can, in principle, be extracted via 
a finite-size analysis of the scaling form of the specific 
heat at low temperatures, which should approach zero 
exponentially with the gap characterizing the decay rate. 
In practice with our QMC method, the specific heat is 
notoriously difficult to measure accurately, especially for 
low temperatures where it is very small. As such, we do 
not present any results on energy gaps in this work, but 
will comment on them briefly below. 

If indeed a 3D photon phase did intervene between 
the two ordered phases, this would be detectable as a 
finite region of phase space wherein both order param- 
eters, along with the energy gap of the system, vanish. 
In Sec. |Vj we will detail precisely which regions of phase 
space we explored in search of this phase. 



IV. N - 7 PHASE DIAGRAM 

In this section, we always take J2 = 0. The goal here 
is to estimate the value of 7 for which the Neel- VBS 
transition takes place for 3 < N < 9. We know from 
previous studies that in the 2D limit (7 = 0), N = 4 
is Neel ordered while N = 5 is VBS ordered; in the ID 
limit (7 = 00), N — 2 is "Neel" ordered while already at 
N = 3, the ground state is VBS ordered. Our main result 
here is that for the 3D isotropic case, in which we expect 
the minimal value of N for which the VBS order onsets 
to be the largest, we find that N — 9 is Neel ordered 
while N = 10 is VBS ordered. Figure [2] summarizes the 
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FIG. 4. (color online). Order parameters for 7 = 1 and 
N = 10 as a function of J2. The squares mark the data 
for the Neel order parameter; circles, data for the VBS order 
parameter. In each case, the data points are obtained at the 
given parameters using a system with L = L z = 16. 



results of a finite size scaling analysis of both magnetic 
and VBS order parameters over the entire N— 7 plane. In 
all cases of 7 tested, we found that the integer values of N 
gave either Neel or VBS order (see Fig. [T] for a depiction 
of the dimer pattern) for the ground state. Our results 
at all values of 7 suggest that the transitions are first- 
order; but without varying J2 (see the next section), we 
cannot be sure since N can only take integer values. The 
results of a detailed study of the order parameters in the 
isotropic case (7 = 1) are shown in Fig. [3] 



V. NEEL-VBS PHASE TRANSITION 

We begin with the isotropic case, 7=1, where we 
have found that the Neel order breaks down for N = 10. 
As shown in Fig. [3j there is clear VBS order for J2 = 
0. With only a small J 2 /J\ ~ 0.028, the Neel order is 
restored. In Fig.|4j we show both of the order parameters 
as a function of J2/Ji- Two conclusions are clear: (1) 
there is no intervening phase between the two ordered 
phases (i.e., no photon phase) and (2) the transition is 
first order. Additionally, our best efforts to extract the 
energy gap from measurements of the specific heat for a 
series of low temperatures suggest that there is indeed a 
gap on the VBS side of the transition for N — 10 and 
J2 = 0, but due to the difficulty of accurately measuring 
the specific heat, as mentioned earlier, we did not pursue 
this study for finite J2 or other values of the anisotropy 
parameter 7 and will not present any results on energy 
gaps here. 

The isotropic case is, of course, very far from the sce- 
nario described in the introduction of weakly coupled 
planes, each in close proximity to where one expects the 
deconfined quantum critical point in 2D. Therefore, we 
now consider 7 = 0.1 where N = 5 is the minimal value 
for which we see VBS order, the same as in the 2D case, 
and this is the value we will study. We again vary J 2 , 
resulting in Fig. [5] Again, the same two conclusions are 
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FIG. 5. (color online). Order parameters for 7 = 0.1 and 
N = 5 as a function of J2 / Ji . Two system sizes are displayed 
to show the trend toward a first-order transition. The smaller 
system size (16 x 16 x 4) is shown with dotted lines and squares 
for the Neel order parameter and circles for the VBS order 
parameter. The solid lines correspond to the larger system 
size with L = 24 and L z = 6. Triangles mark the data for 
the Neel order parameter; diamonds, data for the VBS order 
parameter. The inset shows the evolution of the Neel order 
parameter in Monte Carlo steps (each step shown actually 
corresponds to a bin of 500 steps) at J2/J1 = 0.0245. The 
sharp jumps between values that correspond to each of the 
ordered phases provide additional evidence for the first-order 
nature of the transition. 



clear: no intervening phase and a first-order transition, 
occurring near J2/J1 ~ 0.025. 

There is one subtlety involved with studying the 
anisotropic system, which we shall now discuss. One 
could argue that a sensible system geometry in three 
dimensions is simply a cube, and this certainly makes 
sense in the isotropic limit. However, when the cou- 
pling between the planes, i.e., in the z direction, is 
weakened (strengthened) the system becomes effectively 
longer (shorter) in that direction. In the extreme cases 
where one is very near either the 2D or ID limits, using 
a cubic system actually wastes a large amount of com- 
puter time and obtains data for a point that may be very 
far from the thermodynamic limit. For example, using a 
16 x 16 x 16 system with 7 = 1 is a rather large system 
size perfectly capable of giving insight into the thermody- 



namic limit; however, at 7 = 0.1, this is actually a system 
of 16 weakly coupled planes, each with only 256 sites, far 
less than a typical large system size in a 2D QMC study. 
The situation is even worse for, say, 7 = 20, where one 
has 256 weakly coupled chains of a paltry length of 16, 
which has virtually no chance of capturing large L behav- 
ior. To remedy this problem, we impose different aspect 
ratios of Lj ' L z depending on the value of 7.^ For the 
value of 7 = 0.1, we use L/L z = 4. 

VI. CONCLUSIONS 

We have presented results on the SU(iV) antiferromag- 
net on the cubic lattice. We have connected the results 
of the well-studied 2D and ID cases to the 3D case and 
shown that there is a direct, first-order Neel to VBS tran- 
sition occurring between N = 9 and TV = 10. Addition- 
ally, we have studied a system of weakly coupled planes, 
each in close proximity to their deconfined critical points, 
and shown that while such a system is certainly a viable 
candidate for harboring the much sought after 3D pho- 
ton phase, no evidence of the presence of such a phase 
exists; instead, there is again a direct, first-order transi- 
tion between the two ordered phases. 

Although our exploration of the phase space for the 
model given in Eq. (JlJ is not exhaustive, we have inves- 
tigated the most likely region where one might expect to 
find a stable 3D photon phase. Furthermore, we have 
presented compelling evidence that the continuous Neel- 
VBS transition present in two dimensions immediately 
gives way to a first-order transition upon weak coupling 
of the planes. 

It would be interesting to find other ways to couple lay- 
ers of two-dimensional quantum critical layers that could 
potentially result in interesting three-dimensional spin- 
liquid phases. The identification of such phases in sign 
problem free Hamiltonians will not only allow a detailed 
study of the spin-liquid phase but also permit the unbi- 
ased study of the exotic quantum critical points adjacent 
to such phases. 
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